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I. INTRODUCTION 

Despite the thousands of papers written on the topic 
of spin glasses over the last four decades, even their order 
parameter is still a matter of controversy. The two main 
rival theories are the replica-symmetry breaking (RSB) 
theory of Parisi [l[ , which is motivated by his exact solu- 
tion of the Sherrington-Kirkpatrick (SK) model, and the 
droplet /scaling theory [2]-|4| . 

Much effort has been devoted to the behavior of spin 
glasses in a field in an attempt to discriminate between 
the two theories. According to the RSB picture, there 
will still be a transition in an applied field h, occurring at 
a temperature T c (h) , which defines the Almeida-Thouless 
(AT) transition @. Within the SK model, the AT line 
can be directly calculated. On the low-temperature side 
of the line there is a phase with broken replica symmetry, 
while on the high-temperature or high-field side there 
is a replica symmetric paramagnetic state. Within the 
droplet /scaling theory on the other hand there is no AT 
line: there is no phase transition in any applied field, just 
as in a ferromagnet where the addition of a field removes 
the phase transition. In this theory the low-temperature 
phase in zero field is replica symmetric. 

There have been arguments advanced that the RSB 
picture works for dimensions d > 6 and that the 
droplet /scaling picture only applies when d < 6 . On 
this scenario there should be no AT line for d < 6. Nu- 
merical simulations to investigate this predicted change 
of behavior at dimension six are not feasible. Instead 
investigations have focussed on whether there is an AT 
line in three or four dimensional Ising spin glasses [t| [l(| , 
with no universally accepted outcome. 

In this paper a new calculational technique is used (at 
least it is new in the area of spin glasses) : the 1 /m ex- 
pansion method. With it a strong argument can be given 
that 6 is indeed the dimension at which the AT line dis- 
appears. In particular I argue that the lower critical di- 
mension of the AT line is likely to be the same as the 
lower critical dimension for the spin glass phase of the 
large-m model in zero field. From numerical studies that 
is known to be around 6 [IJ [H| • 



In the 1 /m expansion method one studies first the large 
m, i.e. the m — oo limit, of the m-component vector spin 
glass (m = 3 corresponds to the Heisenberg spin glass, 
m = 1 is the Ising spin glass) and then proceeds to the 
construction of the 1/m-expansion. This procedure was 
advocated some time ago [IJ. Unfortunately, unlike in 
the case of ferromagnets, the large-m limit is not an- 
alytically tractable, and so carrying out the expansion 
requires having a good numerical procedure. While nu- 
merical work for the large-m limit is much easier than for 
finite values of m [14h!1oT ] , it is still non-trivial and this is 
probably the reason why this paper is effectively the first 
paper on the 1/m expansion for spin glasses, although a 
start was made in Ref. (ITsj . 

The large-m, m = oo, model itself is rather unusual. 
Its upper critical dimension is eight [17j ] . That of finite m 
spin glasses is six [l8| . Determination of its lower critical 
dimension has to date not been achieved by analytical 
studies. In her first study Viana [lj| suggested it was 8, 
but a subsequent study revised that to 14 [20]. The con- 
clusion we would draw is that further studies are needed 
to obtain the lower critical dimension analytically. We 
think it is 6, but we would acknowledge that this has yet 
to be properly established. 

What aroused my interest in working on the 1/m- 
expansion was a result of Sharma and Young [2l[ for 
the m-component SK model in the presence of an m- 
component random field. They showed that at T = 0, 
the AT line hits the /i-axis at a finite value, Hat, where 



h 2 



1 



m — 2 



J > 



(1) 



provided m > 2. Thus h AT will have a well-defined ex- 
pansion in 1 jm whose first terms are 



h 2 



J 2 . 



(2) 
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It occurred to me that it might be possible to obtain 
an expression for h 2 AT outside the SK limit in terms of 
correlation functions of the large-m model, by means of 
a 1/m expansion. 

In the large-m limit itself there is no AT line, i.e. there 
is no transition in the presence of a m-component random 
field, even in the SK limit. There is a transition in zero 
field and the low-temperature phase is replica symmetric 
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[13|. In order to get an AT line one has to go to order 
1/m. We show in this paper that Hat at T = is given 
by 



1 



1 



(3) 



where xu is the local susceptibility at site i in the large- 
m limit. In the SK limit, xu — 1/*^) so our general 
expression in Eq. ([3]), which should be valid in any di- 
mension where there is an AT line, is consistent with that 
of Sharma and Young for the SK model. 

Sharma and Young [2l| found also that the degeneracy 
of the fields which go soft at the AT line in the replica 
field theory of the m-component field spin glass model 
is n(n — 3)/2, for all values of m, just as for the Ising 
spin glass at its AT line .6], This suggests that the AT 
line in the m-vector model is in the same universality 
class as that of the Ising model and that it will have 
the same effective field theory as was used in [6j. Thus 
although our work is focussed on the AT line in the m- 
vector spin glass with m large, our conclusions about its 
lower critical dimension are applicable also to the lower 
critical dimension of the AT line of the Ising spin glass. 

Unfortunately, to use Eq. one needs information 
on the xa °f the large-m model in dimension d > 6, 
which would be very challenging to obtain numerically. 
Fortunately, there are arguments which enable one to 
estimate their magnitude: xu ~ l/?c, where T c is the 
transition temperature of the large-m model in dimension 
d. Then 



T c {df/m. 



(4) 



We suspect that the lower critical dimension of the large- 
m limit to be six, so T c should go to zero as d — > 6. Then 
the AT field h 2 AT should also go to zero in six dimensions: 
the lower critical dimension of the large-m limit and the 
AT line should be the same. 

In Sec. |ll]the model which we shall study is defined and 
the formalism for carrying out the 1/m expansion is set 
up. It is a saddle-point procedure and the leading terms 
for the large m limit are also given in that section. The 
formalism is complicated because of the need to obtain 
the spin glass susceptibility, which is the susceptibility 
which diverges at the AT line. In Sec. IIIII the one loop 
calculation is done for the energy of the system and the 
spin glass susceptibility at zero temperature, and from 
the latter we obtain the remarkably simple formula for 
the AT field at zero temperature given in Eq. Our 
arguments on the form of xu are given in Sec. |IV] We 
conclude with a discussion in Sec. [V] Finally in the Ap- 
pendix, it is explained why the problem of understanding 
the AT line is so hard: it is because it is fundamentally 
a non-perturbative problem. 



II. FORMALISM FOR THE LARGE m LIMIT 

We shall study the m-component spin version of the 
Edwards- Anderson spin glass model in an m-component 
random field. This has Hamiltonian 



u 



^ N m N m 

~ ^ Jij ^ SifiSjft — ^ ^ fli^S; 



(5) 
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where the random field components hi^, [i = 1, 2, • • • ,m, 
are each drawn from a Ga ussi an distribution of zero mean 
and variance h 2 so that /i? = h 2 . The m-component 
vector spins Si are chosen to be of fixed length: 

, 1 = 1,2,..- ,N. (6) 



E 



/ j Sip 



The spins will be taken to sit on the sites of some 
d-dimensional lattice where the sites are labelled i = 
1, 2, • • • ,N. The coupling between the spins, Jij, could 
have any distribution which leads to a spin glass phase. 
The form of this distribution is not crucial to what fol- 
lows but it is convenient to imagine that it is between 
nearest neighbor sites on the lattice and has a symmetric 
Gaussian distribution of width J and zero mean. The 
SK limit is when the interactions are between all pairs of 
spins: the width of the distribution has then to be set to 
J/\/~N in order to make the energy in the thermodynamic 
limit extensive. 

Because we want to study the AT line, we need to 
calculate the replicon (spin glass) susceptibility xr- It 
is the divergence of XR as the field h is reduced while 
keeping the temperature T fixed which determines the 
location of the AT line h(T) [5(. In order to calculate 
the replicon susceptibility it is convenient to employ two 
copies of the system: the Hamiltonian then is 



H[S,T] = 



^ N m 

2 E/ ^ y~^('~V'5j> 

N m 



EE 
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+ Til 



(7) 



where the fixed length spins Tj are again such that 
Tfn = to at all sites i. The Edwards- Anderson order 
parameter q can be expressed in terms of and X^: 

^ V ra _^ A m 



zij,T ifJj }. (8) 



(The overline indicates the average over the random field 
distribution. We shall have no need to average over the 
bond distribution of the couplings Jij). Similarly the 
replicon susceptibility \R is 



1 



Xr 



Nm 



0- 



! E 

*iMi.J> 



1 



[{SinSju) — (Sifj,)(Sj U )} 2 — ]Ym^ 



[{SiM - (Si,i)(S jv )][(T ip T jv ) - (T ip )(Tj v )], 
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where here /? = 1/T. 

To handle the constraints imposed by having fixed 
length spins, it is convenient to make use of the rep- 
resentation fl3l. l22j] 

rn „ioo qjtt m 

S(m - £ Si) = / exp[/3(m - £ 5^)^/2]. 

J —too 



We have introduced the variables h u and /i m in order to 
generate various correlation functions from differentiat- 
ing the free energy, F/T = — \nZ 2n , but after differen- 
tiating with respect to them one sets at the end of the 
calculation 



hD = h M = h. 



(9) 



The partition function of two copies, Z 2 , where Z is 
the partition function of a single copy derived from the 
Hamiltonian of Eq. (JS|), can be written 

--/n^/n^£n*£n^ 

expi -/3H[S,T] + 



The free energy F then of the replicated system is just 
2n times the free energy of the unreplicated system. For 
example, given F one can obtain the Edwards- Anderson 
order parameter q as defined in Eq. (j5J| via 



E 



/3(m - Sl)Hi/2 + /3(m - £ 1^/2 

M— 1 M— 1 



In order to calculate the free energy and correlation 
functions such as XRi averaged over the random field s, it 
is necessary to use the replica trick and obtain Z 2n . As 
usual n has to be set zero at the end of the calculation. 
To this end one attaches replica labels so that Hi, Ki, 

and T v become Hf, Kf, and T? u , where a 
1, 2, ■ • • , n. Then 



n 
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On averaging over the random fields 



ex p { I £ £ J « £( s s. s ?, + TC) 
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Second derivatives with respect to fi 2 h 2 M yield the repli- 
con susceptibility: 



8 2 {F/T) 
d((3 2 h 2 M ) 2 
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= n 2 l^S^) - (S^)(S^}mM - (T V )(T ]U )] 
+0(n 3 ) = n 2 NmxR/(3 2 as n -> 0. (11) 



We thus will need to calculate the free energy F to order 
n 2 to get xr- 

The next step is to integrate out the spin variables 
S"„ and T°j. This is readily done using a cumulant 
expansion. It is convenient to let the replica indices 
on H" run 1,2, ■■■ ,2n, and set Kf = Hf when a — 
n + 1, n + 2, ■ ■ ■ , 2n. We shall define 



At 



and call its inverse 



Y a . = r — 1 1 ■ ■ 



Then the remaining integrals for Z 2n are 



(12) 



(13) 



Z 2r. 



n 
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exp[mS({H?})], (14) 
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where, up to the third cumulant (and terms of order n 3 ) 

a 2n i 
S({H?}) = P - £ \ H i + -s lndet(xg-/j8) + J&Xg] + 



?' a — 1 



27J 



n 2n 



t +^ E +<E E 



a /3 
Xij Xji 



2n 



E +$> E + 

\ a,/3,7— 1 a,/9,7— n+1 

n 2n n 2n \ 

3*2*1, E E +3^>4E E xgx&xz,. 

a,j3=l 7=n+l a=l /3,7=n+l / 

where we have adopted the notation that repeated site 
indices are to be summed from 1 to N. 

Because of the factor of m, which we are taking to be 
large in Eq. (fl4| . the integrals over H" can be done by 
the method of steepest descents. The leading term can 
be found from the solutions of 



dS 
dH? 



0, for i = 1, • • • , N, and a = 1, 



There is a solution of these Nn set of equations which 
has replica symmetry, Hf = Hi for a = 1, • • • , n. The 
values of Hi can be obtained by solving the N saddle- 
point equations, 



P - Xm ~ P h D E XikXki ~ n/3 2 (h% + h%[) ^ XikXklXli 



k.l 



XikXklXlmXr 



(15) 



fc.J,'. 



We shall work in the low-temperature limit T — > or 
j3 — > oo, which permits a number of simplifications. In 
this limit, the term in \a, which is finite at T — 0, can 
be dropped as it is negligible in comparison to the terms 
with factors of (3. In the limit of n — > 0, Hi — > H® which 
can be determined by solving the N equations 

1 = h D^2xikXki- 

k 

Using the saddle-point solution the free energy to order 
n 2 in the large m limit as T — > is then 



F/T = -pnm ^ Hi 



(3 2 n 2 mN 



2h 2 D 



K + h A M ). (16) 



It follow on using Eq. (TT0|) that q = h 2 M /h 2 D which equals 
1 on setting hp = h = Km- This is what would be 
expected at zero temperature. The replicon susceptibility 
can be calculated using Eq. (fTTj) and is given by \R = 
1/hl = 1/h 2 . 

It will be important for what follows at one-loop order 
that we know the leading changes which a finite value of 



n would make to the Hi. If we set H k -> H { k 0) + AH k , 
then 



X^X^-E^AH.xg- 



(17) 



From Eq. (|15[) . to order n, 

Ex* = rr ~ P n ( h D + hit)£r, 



where 



A« ~ E XikXklXli ^ ■ 



(18) 



k.l 



By substituting for \ik using Eq. (1171) one has to order 



^X^^H hX tlxf=-Pn[h% + K 
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m)tt- 



which has solution 



AH k = A H 



pn 
21% 



(h^ + hlf), k = l,---,N. (19) 



S({Hi}) goes like n as n —> 0. Hence changes of order n 
in Hi will only contribute terms of order n 3 to S. This is 
because 5 is stationary with respect to variations of Hi . 

We shall need a special case of Eq. (fT7|) in the one- loop 
order calculation: 
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III. ONE-LOOP ORDER CONTRIBUTION 

The one-loop order contribution to F comes from the 
Gaussian fluctuations about the saddle-point solution 
H? = H h Writing 8Hf = H? - H u the terms in the 
expansion of S in Eq. (fT4|) to quadratic order in SHf are 

2n 



111 E 5H?M aP (ij)SHf, 



ij a,/3— 1 

where the matrix 



d 2 s 



M « m = BH?mr 

t J 

is evaluated at the saddle-point Hf = Hi. 

M a p(ij) = SapCij 

+(3 2 h%D tJ + 2n^{h% + f&hfaPij, l<0<n, 
+[3 2 h 4 M D ij + Anf3*h 2 D h A M P l3 , n + 1 < p < 2n, 
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provided a < n. The case of a > n + 1 is obtained by 
switching around the two forms dependent on 8. Here 



D, 



Si 



cy 



^ ^ XikXilXkjXlj i P%j ^ ^ X«fcX/c£X/jXjmXrmi 

k,l,rn 

^ ^ XijXjkXklXli-! Qij ^ ^ X«jXjfcX/c/X/mXmi- 

k,l.rn 



k 

2n 2 B 3 (h 6 D + 3h 2 D h i M )Q l 



The matrix M a p(ij) is a 2n x 2n matrix in replica space 
with Cij along the diagonal. In the two n x n blocks near 
the diagonal the elements are all B 2 h' D Dij + 2B 3 (h D + 
h 2 D h\ { )Pij, while in the two nx n blocks off the diagonal 
the elements are all B 2 h\jDij + A^fijjh^Pij . Such a 
matrix has a 2n — 2 fold degenerate eigenvalue equal to 
and two non-degenerate eigenvalues 



Uij Cy 

Vij = Cij 



nB 2 (h D 

d- 



h 4 M )D. 



nB 2 (h 4 
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D n M)Piji 
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If we denote the one-loop contribution to Z 2n by exp S± , 
1 



51 = -| [(2n- 2) In det Cy 



In det Uij + In det Vij] . 



To make analytical progress we shall work in the limit 
T — > 0. Furthermore we will find that if we keep h 2 D and 
h 2 M small, of 0(J 2 /m), the expected magnitude of the 
field at the AT line, then substantial simplifications are 
possible. 

Let us begin by simplifying Ca with the aid of Eq. 
(TTSl). It becomes 



C l% = 2(/3 - XiijXii, 

if one uses the identities Su = xu J2k,i XikXkiXki&nd 
Qu = Xii2~2k,l, m XikXklXlmXmi- As T -> 0, Xu remains 
finite, so C u -t 2Bxa- 

The off-diagonal terms in Cy with i ^ j and n = 
are in the T — > limit 2Bh 2 D Xij Ylk XikXkj- Thus in the 
limit we shall work, when h 2 D ~ J 2 /m they are negligible 
compared to the diagonal terms and we can write Cij — 
2j3xu&ij- With this approximation, it is straightforward 
to evaluate 51 to 0(n 2 ). 
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nBh A D ^2 Dn/Xt 
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For the terms linear in n we have to use for Da its form 
correct to order n to get the n 2 terms in 51: 



Di 




1 



[1 - 2nB(h 



and also use Eq. (|20l) for xu- Now 

Pit = ^2 Xim XikXklXU = TT A *> 



k.i 



and the double sum 



as for small hjj it is dominated by the term with i — j. 
The term in 51 equal to — n 2~2i ln(2/3x^) is subdominant 
as T — >• 0. Assembling all the terms together the leading 
terms up to one-loop order are 

FIT = -Bnm V Hi + -Bn V — 

V ^— ' V.-.- 



B 2 n 2 mN 

2h 2 D 
B 2 n 2 {h\ 



(h% + h 4 M ) 
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2h D h M h\f) 



8h% 



x 2 - 



(21) 



In this equation, xu an d -Hi are &t their n — ?• values. 

By using Eq. (fTQ|) one can confirm that this expression 
for F/T still implies that q = 1. By using Eq. (TTT]) the 
replicon susceptibility is, after setting ho = h = h,M, 



1 11-^1 

XR ~ W + mh i N^xJi 



(22) 



Terms which have been neglected, like the terms of 
with i j^z j give rise to subdominant corrections of order 
l/(mh 2 ) in xr which are negligible compared to the sec- 
ond term in Eq. (I22[) provided we are only interested in 
values of h 2 ~ J 2 /m. 

In the SK limit of the large m model, Xu h as n0 site 
dependence in the limit N — > oo and at T = equals 
1/(J 2 + h 2 ) 1 ! 2 . If we take h? to be small, say, of order 
J 2 /m, then for the SK model Eq. (|22|) reduces to 



XR 



1 



J 



1 



ih 4 



(23) 



with h 2 AT 



J 2 /m. This result is therefore in perfect 
agreement with that of Sharma and Young [2l[ for the 
value of the critical field at T = 0. For all d > 6 one is in 
the classical regime where the exponent 7 characterizing 
the divergence of \R as one approaches the AT line is 
equal to 1, a feature which is captured by the simple 
pole expression of Eq. (|23|) . 

Outside the SK limit, the expression for the AT field 
at T = generalizes to 



h 2 



=—y- 



(24) 



We shall discuss the consequences of this for the proper- 
ties of spin glasses in Sec. IIVI 
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The ground state energy E of the system is obtained 
from Eq. (|2ip by dividing the linear term in n of F by 
the total number of replicas 2n, so 



In the SK limit, Hi is site independent and 
J 2 



(25) 



Hi = 



V J 2 + h 2 



+ ^J 2 + h 2 . 



(Solution of the m — oo SK spin glass model is aided by 
the observation of de Almeida et al. [23J that it is replica 
symmetric and has the same free energy as the spherical 
spin glass, [24], although the physics of the two are quite 
different |13|). 

The 1/m expansion done in this paper will only have 
utility in the region of fields and temperatures above the 
AT line. Nevertheless, provided that we take h 2 > h 2 AT , 
but still of order J 2 /m, it is possible to obtain the ground 
state energy at h = correct to order 1/m, as can be seen 
by expanding out the expression for Hi in a power series 
in h 2 : H = 2 J + h 4 /4 Hence in the SK limit, the 
ground state energy is given by 



E/m 



-J 



J 

4m 



0(J/m 2 



(26) 



Setting m = 1, the ground state energy per spin of the 
Ising spin glass is estimated to be — 0.75J. It is actu- 
ally equal to —0.763... J, [lj], suggesting that the 1/m 
expansion might have utility even for Ising systems! 

In general one would expect that the 1/m expansion 
would only work when the state about which one is ex- 
panding, the large-m limit, is similar to the state at finite 
values of m. The low-temperature phase of the SK model 
in a random field at finite values of m has full replica 
symmetry breaking and is quite different to the param- 
agnetic, replica symmetric state of the large-m limit. So 
it is rather surprising that putting m = 1 gives quite a 
good result for the ground state energy of the SK model. 
However, in say two dimensions, there is certainly no AT 
line at finite values of m, and one then could hope that 
the 1/m expansion for the energy should work well. It 
should be easy to obtain by numerical methods the Hi 
and the xu m both two and three dimensions which have 
to be inserted into Eq. ([231) in order to check this out. 



IV. THE AT FIELD IN FINITE DIMENSIONS 

In order to calculate h 2 AT using Eq. (f2"4"|) one needs 
information on the Xiii which are, alas, unknown. We can 
only make progress by examining their likely dependence 
on the value of T c . 

So long as there is a finite temperature transition in 
the large-m limit model, \a wm be finite as h goes to 
zero and it is its value at h — and at T = which 



should be used in Eq. ([24"|. At h 
equations simplify to f3 = xu [H 

T<T r , 



- 0, the saddle point 
when T > T r . For 



Xit = 13(1 - Qi), 



(27) 



where 



1 m 

711 < J 



Below T c , there is a Bose-Einstein like condensation into 
a large number of modes, mo. For the SK limit mo ~ 
N 2 l b [13]. Each component /i of (Sip) is proportional to 
a different one of the mo nearly null eigenvectors of the 
Aij matrix, qi is the sum of mo terms 



m 

m f-~* 



The site average of qi is the Edwards- Anderson order 
parameter q. 

In the SK limit q = 1 - T/T c , so x = P(l - q) = 1/T C 
for all T < T c . Outside the SK limit, q is not known 
precisely, but as T — > it would seem natural to expect 
that it goes like q = \—aT/T c , where a is a constant. This 
is certainly consistent with the results in Ref. [l6| (albeit 
on the long-range one-dimensional long-range version of 
the large m model) . Similarly one would expect that the 
low-temperature form of the qi should be qi — l — aiT/T c . 
Then Xii = a i /T c at T — and 



h 2 



AT 



m T2 

™ c ' 



(28) 



Not much is known about the coefficients a^. In the SK 
limit, Xu = 1/J = 1/T C so a t = 1 at all sites. Thus 
in that limit k = 1. We shall assume that k(d) remains 
finite as d — > 6. 

Normally, as a system approaches its lower critical di- 
mension its transition temperature falls to zero. As we 
suspect that the lower critical dimension of the large-m 
limit is 6, it would follow then from Eq. (|28|) that the 
AT line should also disappear as d — > 6. 

The dependence of T c on dimensionality d can be elu- 
cidated further if the McMillan RG equation Q] applies 
to the large-m model. Certainly it seems to work in 
the one-dimensional version of the model studied in Ref. 
16]. McMillan argued that the flow of the temperature 
T under a renormalization group transformation near the 
lower critical dimension where 9 is small was as 



dT 
d\nL 



= -9T + cT 3 /J 2 



(29) 



and c is a constant of order one. Provided that 8 is 
positive, this has a fixed point when T 2 — 6 J 2 /c. Near 
the lower critical dimension 9 changes sign, so let us write 
it as 9 ~ (d — 6). Then close to six dimensions, 



h 2 



d-6 



m 



k(d)J 2 



(30) 
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Thus according to this calculation, at one-loop order 
the AT line should not exist below 6 dimensions, or more 
conservatively, the lower critical dimension of the large-m 
model. 



V. DISCUSSION 

The remaining major task is to show conclusively that 
T c , the transition temperature in the large-m limit at zero 
field, does indeed go to zero as d — s- 6 or in other words, 
that six is the lower critical dimension in the large-m 
limit. There are intriguing features of the large-m limit 
which suggest that this problem is not as hard as that of 
directly determining the lower critical dimension for the 
AT transition. There is a suggestion in [I?) that the ex- 
ponents in the large-m limit can be related to the critical 
exponents of a problem without disorder in two fewer di- 
mensions, which is an example of dimensional reduction. 
Dimensional reductions can often be related to some su- 
persymmetry, so maybe an elegant demonstration that 
six is the lower critical dimension will be possible. 

Assumptions have had to be been made in this cal- 
culation because the xu are n °t known outside the SK 
limit. We could make some progress by using the usual 
assumptions as to the form of q at low temperatures, but 
we could not work out the xu explicitly. Unfortunately, 
numerical studies of the large-m model near six dimen- 
sions to check the assumed form of xu are unlikely to be 
practical. 

However, it is possible to do numerical work on the 
one-dimensional version of the model with long-range in- 
teractions, as in Ref. [l6j]. Here the interactions between 
the spins are still random but fall off with a power of their 
separation ry: Jy ~ e ij/ r fj- The are Gaussian ran- 
dom variables of variance J 2 , while the exponent a can 
be tuned to mimic various dimensions of the short-range 
model in d dimensions [IB]. In zero field, the SK model 
corresponds to a < 1/2, the upper critical dimension 8 
of the large-m limit maps to a = 5/8. The lower critical 
dimension (which I think is 6) corresponds to a = 3/4. 
The numerical work in Ref. 16] provides good evidence 
that T 2 ~ (0.75 — a) J 2 . 

Now I believe that the AT line at any finite value of 
m disappears for d < 6 or values of a > 2/3 7]. But 
our 1/m expansion to first order in 1/m would predict 
that h 2 AT is finite until a reaches 3/4. I suspect that 
non-perturbative effects will remove the AT line in the 
interval 2/3 < a < 3/4. This will be discussed in a 
separate paper on the AT line for large but finite m in 
the one-dimensional model [25j ] . 

It is easy to see that non-perturbative effects must be 
important. Suppose there were an AT line for 2/3 < a < 
3/4. There is no stable perturbative fixed point in the 
vicinity of the "upper critical dimension" , i.e. when a > 
2/3, p}, indicating that the calculations in this interval 
have to be non-perturbative. Non-perturbative effects 
must also be important for a < 2/3 in determining the 



value of h\ T , even though the critical behavior at the 
AT line is there controlled by the Gaussian fixed point 
0. Non-perturbative effects must reduce h\ T to zero, 
probably as some power of (2/3 — a), so that h 2 AT joins 
smoothly to being zero for a > 2/3. 

For the spin glass model in d dimensions it is a happy 
coincidence that the lower critical dimension of the AT 
line and that of the large-m limit both seem to be six. It 
is likely though that the same type of non-perturbative 
effects which must dominate the one-dimensional model 
for a close to 2/3 might affect some of our estimates of 
h 2 AT when d is close to 6. Our one- loop order calculation 
has suggested that h 2 AT ~ T 2 A(m,d) where A(m,d) ~ 
1/m as m — > oo. Contributions from higher orders in 
the loop expansion would still be expected to go as T 2 , 
and would just add further terms in the 1 /m expansion of 
A(m, d), as in Eq. @ for the SK limit (recall that in this 
limit J 2 = T 2 ). Some of the dimensionality dependence 
of h 2 AT will arise from that of T 2 , which I think varies 
as [d — 6) J 2 . Non-perturbative effects might at fixed m 
make A{m, d) vanish as some power of (d— 6) just as they 
might produce powers of (2/3 — a) in the one-dimensional 
long-range model. But even without being able to solve 
the non-perturbative problem, we can still be confident 
that there is no AT line in less than six dimensions, as 
T c will vanish in six dimensions. 

Finally, the 1/m expansion procedure is quite general 
and can be used for any quantity in principle. It would 
be good to apply it to problems for which there is no 
analytical work, such as the value of the droplet exponent 
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Appendix: Why calculating the AT line near six 
dimensions is hard 



In this Appendix we explain why investigating the be- 
havior of the AT line as d — > 6 is hard for Ising spin 
glasses. It is because it is a non-perturbative problem. 

The Ginzburg-Landau- Wilson free-energy functional 
for the Ising spin glass written in terms of the replica 
order parameter field is TSj 



F[{Q a p}\ = fd d x [|r£ a</3 Q^ + §£ a</J (VQ o/) ) 

+ IT Ea<)5<7 Qa^Q^Qja ~ h 2 J2a</3 Qa/3 + 0(Q 4 ) 



where h is the applied field. A simple scaling analysis of 
the terms in the functional shows that the natural size 
of h 2 AT is ~ |r| 2 /u>; this remains the correct scaling form 
for all d > 6. Perturbation theory is in terms of the 
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the limit when x is large or at least fixed. This is when 
T is fixed (at a value less than T c ) or \r\ is set at some 
large value. But this is the regime which is basically non- 
perturbative and no calculations valid for it have been 
obtained: the only calculations which one can be done 
systematically are those which are perturbative. It is 
this fact which makes determining the form of the AT 
line in the limit d — > 6 so challenging. 

Furthermore, the general problem of understanding 
how one can have a replica symmetric state for d < 6 
is also a non-perturbative task |8|, which explains why 
analytical progress has been so slow. 
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